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ABSTRACT 

Boundary layer velocity and temperature profiles were 
measured in nitrogen near its thermodynamic critical 
point ftowing turbulently upward over a heated flat plate 
with constant surface heat flux. The measurements were 
made with a combination thermocouple and pitot-static 
I probe. 

Data were taken with the fluid in the boundary layer 
spanning the critical temperature and also at the same op- 
erating conditions except with the boundary layer at tem- 
^ peratures away from the critical temperature. These re- 
sults showed that the near-critical vebcity profiles were 
> quite different, exhibiting a maximum between the wall 

( and free stream. The near-critical profile had a much 
larger Grashof number suggesting that a strong body force 
effect was present. The data are examined in terms of 
the velocity and temperature profile response to variations 
in the system parameters. Proximity to the critical point 
appears to exaggerate body forces even at quite high 
Reynolds numbers. 

I NOTATION 

I D tube diameter, mm 

g acceleration of gravity, cm/sec^ 

Gtx Grashof number, p^(p^ -p^)x3g/^^ 

h heat transfer coefficient, W/cm^ °K 

L tube length, mm 

P pressure, atm 

q heat flux, W/cm^ 


r} dynamic viscosity, gm/cm-sec 
p density, gm/cm^ 


RCj^ Reynolds number, 


T temperature, 
u velocity, m/sec 

X distance from leading edge of plate, cm 
y normal distance from the plate, mm 


Subscripts 


b fluid bulk conditions 
c critical point conditions 
w wall conditions 

X conditions at axial position x along the plate 

00 free stream conditions 

Superscript 

^ transposed critical point conditions - temperature as a 
function of pressure at which the constant pressure 
specif heat curve maximizes 

INTROPUCTION 

Because of the strong, nonlinear thermal property 
variations of fluids near the thermodynamic critical point, 
heat transfer, and fluid flow of near-critical fluids are 
challenging research an^ design problems. In recent years 
the necessity of oper*>* ag near the critical point has com- 
plicated the design analyses of hydrogen rocket engines, 
high pressure steam generating systems and or 3 romagnets. 
In the future, the possible use of supercritical methane as 
an aviation fuel may present a similar problem. 

The field has been actively researched for about 
25 years. This research has been comprehensively sur- 
veyed in recent articles by Hendricks et al. (1), Hall et al . 
(2), and Petukov (3). It is apparent from these surveys 
that there exists a lack of data of sufficient detail to make 
an accurate description of the velocity and temperature 
fields, especially in forced convection. This in has 
hampered the selection of a realistic model for analysis. 
The literature does not lack suggested analyses. What is 
lacking are the velocity and temperature profile data needed 
to definitively test a given flow model. 

Only two experiments have been performed in which 
profile measurements have been made during heat transfer 
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to a near-critical fluid. Wood (4, 5) and Wilaon (6) both 
made velocity and temperature mirveys near the exit of a 
smooth tube with a constant heat flux boundary condition. 
Wood (4» 5) flowed carbon dioxide vertically upward in a 
22.9 mm diameter tube. He made measurements at L/D 
B 30.7, and could survey to within 0.25 mm of the wall. 

He observed that under conditions of T|^ < T’*' < the 
temperature profiles flattened considerably and heat trans- 
fer coefficients rose sharply as the bulk fluid approached 
the critical temperature. He also observed that the velo- 
city profile achieved its maximum value between the wall 
and centerline, yielding a characteristic M-shape. His 
experiment was hampered somewhat by the fact that the 
velocity and .emperature profile measurements were made 
separately, requiring matching of two runs. At present 
no analyses of Wood's data have been published. Wilson (6) 
flowed hydrogen, also vertically upward, at L/D = 128.6 
in a 35.0 mm diameter tube. Wilson could not make 
measurements as close to the wall as Wood; however, he 
employed a combined thermocouple, hot wire, and pitot 
probe and took continuous data. Wilson's results con- 
firmed the presence of the M-shaped profile and exhibited 
a great deal of noise in the profile traces. Wilson's re- 
sults also showed a strong spike in the waU temperature 
distribution axially along the tube (see also (1)). Wilson 
correlated his data using a pseudo -boiling model. 

The potential of using the hot wire in the near-critical 
region has also been examined by Bourke et al . (7) in adi- 
abatic tests with C02- Despite our lack of knowledge of 
the heat transfer processes in this region the hot wire 
shows real promise for meaningful measurements. 

In addition to these velocity and temperature survey 
experiments, Sabersky and Hauptmann (8) published some 
outstanding color photographs of a near-critical carbon di- 
oxide boundary layer on a horizontal flat plate. These 
pictures coupled with the data, including hot wire signals 
measured at 0.25 mm above the plate, gave some insight 
into the nature of the fbw. They concluded that the flow 
could be treated as single phase turbulent ^ow, possibly 
requiring modifications for enhanced turbulent fluctuations. 

Recently several investigators (9-14) have directed 
their attention to the phenomenon of the large increase or 
"spike" that occurs in the wall temperature under con- 
stant surface heat flux conditions in near-critical fluids 
Rowing in smooth tubes. It was found that the "spike" 
could be attenuated mechanically (e.g. , swirl devices) (10, 
13, 14) or by the presence of a strong body force effect 
(9-13). While the wall temperature spike is not studied in 
foe present experiment, strong body force effects did oc- 
cur and the above investigations are related accordingly. 

In these tube experiments (9-13) foe body force effect was 
demonstrated by flowing both upward and downward. It 
was determined that body forces would influence the wall 
temperature if foe Grashof number divided by foe square 
of the Reynolds number were greater than 10"^. The nu- 
meric value will of course depend on the definition of 
Reynolds and Grashof numbers; it is the parametric 
grouping tha; is o! interest. 

Many analyses and flow models have been directed 
toward ne«r-crltl<*al heat transfer. The present results 
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should bo of particular interest in testing the various modi- 
fications of mixing length theory proposed over vhe years 
(15-21). Of special interest is the work of Hsu and Smith 
(17) because it included the body force effect and predicted 
the Grashof-Reynolds relationship mentioned above. 

The present study was undertaken to obtain velocity 
and temperature boundary layer profiles in a geometry 
which was as simple as possible. This was done by flowing 
Tiear-critica" nitrogen vertically upward over a flat plate 
having a constant heat flux at the surface. The effort here 
was to avoid some of the complex phenomena that occur in 
tubes - particularly the wall temperature spike. Wilson (6), 
for instance, reported the wall temperature spike, whereas 
Wood (4, 5) did not. U was felt that the simpler geometry 
would be more amenable to analysis. Furthermore, the 
flat pl'^te configuration allows more parametric control. 

For instance, the fluid free stream properties can be con- 
trolled independent of the surface heat flux, which is not 
the case in tubes. It should be emphasized that the tube 
data are very important and closer to reality. The present 
work intends to complement foat of Wood (4, 5) and Wilson 
( 6 ). 

The range of operating conditions obtained are indi- 
cated on Fig. 1, the density curve for near-critical nitto- 
gen at P/Pq = 1.10. The thermal properties were obtained 
from a computer program of Pew et al , (22), which em- 
ployed the equation of state of Strobrid^e (23) and the vis- 
cosity correlatfon of Brebach and Thodos (24). Of all foe 
strongly varying properties associated with foe critical 
point, it is the density which appears to play foe most im- 
portant role. In the present e^qiOriment data were acquired 
at operating conditions such that: (1) bofo the waU and free 
stream were below the critical temperature, (2) bofo the 
waU and free stream were above the critical temperature, 
and finally, (3) foe waU was above and the free stream was 
below foe critical temperature. The data presented herein 
are bounded by the following ranges of foe parameters. Of 
course, not all combinations were covered. 

PyPc«1.10 

0.70 < T„/T„ < 0.99 
00 c 

0.80 <T^/Tc <1.75 
0 < q < 15 W/cm^ 

7 . 0x10^ < Re^^ < 17 . 6X10® 

Pq b 33.5 atm, Tq b 126.2^ K, » 0.311 gm/cm? 

The data will be presented in a manner as to illustrate 
the influence of parameters such as: free stream tempera- 
ture, surface heat flux, Reynolds number and Grashof num- 
ber on foe boundary layer profiles. 

EXPERIMENTAL APPARATUS 

The experiment can best be described by thinking of it 
in terms of a conventional wind tunnel, except here the 
fluid is near-critical nitrogen flowing at about 1 meter per 
second. The test apparatus is illustrated schematically in 
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Fig, 2, It was a onee-through or blow-down type system. 
Liquid nitrogen was transferred from a low pressure sup- 
ply to the 0.4 high pressure dewar, The fluid temper- 
ature was controlled by bubbling nitrogen gas through the 
dewar. Nitrogen gas, at the desired pressure, was applied 
to the top of the fluid to force the fluid through the system. 
System pressure and flow rate were maintained by the use 
of throttling valves. The system was heavily insulated and 
was precooled from an auxiliary supply. The once-through 
system avoided the problem of pump oscillations, but not 
the system oscillations characteristic of near-critical 
fluids (1). 

The test section was of rectangular cross-section of 
about 6 to 5 aspect ratio. It contained two 2. 5 cm wide by 
15 cm long nichrome plates back to back. The first centi- 
meter was a copper electrode and was unheated. The pri- 
mary surface was instrumented longitudinally with 11 
chromel-constantan thermocouples electrically isolated 
from the heater. The opposite surface served as a guard 
heater and was also instrumented. Heat was generated by 
the dissipation of direct- current power. Figure 3 is a 
photograph of the heater plate assembly. 

The boundary layer surveys were made with a combi- 
nation pitot- static and thermocouple probe which was tra- 
versed through the boundary layer toward the primary 
heater surface. The distance from the primary heater 
surface to the opposite wall was about 2.4 cm, which was 
approximately 10 times the boundary layer thickness com-» . 
puted from a lAth power profile. On the back side the 
distance to the wall was 0.72 cm, about 3 times the com- 
puted boundary layer thickness. The temperature portion 
of the probe was an exposed ball, chromel-constantan 
thermocouple encased In a sheath with an overall diameter 
of 0. 20 mm. The dynamic pressure part of the probe was 
a two leg pitot-static probe of rectangular cross-sdction at 
the tip arid had an overall tip height of 0. 20 mm. Pictures 
of the actual probe are shown in Fig. 4. The total pres- 
sure leg had a 0.075x0. 70 mm opening at the tip. The sta- 
tic port was a 0. 25 mm diameter hole located six thick- 
nesses from the tip to avoid separation. Measurements 
could be made to within 0.10 mm of the wall. All of the 
surveys reported herein were made 13. 1 cm from the lead- 
ing edge ef the plate. 

An important feature of this design was making both 
legs identical. This was done in order to match the re- 
sponse lags in the lines, reducing the noise, and also to 
assure the same density gradient in each leg to avoid sta- 
tic head errors. The probe noise was assessed to be due 
to very small static pressure perturbations. This is dis- 
cussed in detail in (25). Wood (4,5) and Wilson (6) both 
observed probe noise. Wilson's level of noise was about 
the same order as in the present experiment. 

Flow rates were measured with venturi flow meters. 
All system pressures were measured with strain gage 
transducers, except at the probing station where the sys- 
tem static pressure was measured with a high precision 
transducer accurate to fcO.0.5%. The dynamic head was 
measured with a variable reluctance type transducer. 
System temperatures were measured with platinum re- 
sistance thermometers. The heater plate thermocouples 
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were staniard chromel-.x>nstantan. All signals, including . 
test section voltage and current, were sensed with a pre- 
cision digital voltmeter and recorded on a digital data ac- 
quisition system. All measurements were checked by some 
kind of redundant instrument. The accuracies of the major 
parameters based on percent of the average value measured 
are presented in Tabic I. 


Table I. - Accuracy of Major Parameters Based on the 
Average Value Measured 


Free stream temperature 

±0.2% 

±0.2® K 

Free stream pressure 

±0.1% 

±0.04 atm 

Free stream density 

±1.0% 

±0.006 gm/cm^ 

Free stream velocity 

±1.6% 

±0.02 m/seo 

Mass flow rate 

±1.3% 

±0.008 kg/sec 

Heat flux 

±1.0% 

±0.1 W/cm^ 

Wall temperature 

±1.0% 

±2® K 

Probe temperature 

±0.6% 

±0.6® K 

Dynamic pressure 

±2.0% 

±0.001 psi 


The manner in which the data were taken has a 
bearing on the results, particularly the number of points 
which could be acquired during a probe traverse. The Sys- 
tem was brought up to nominal operating conditions and al- 
lowed to stabilize on some near-nominal operating condi- 
tion. At this point a traverse of the probe from the free 
stream to the plate was begun. The traverse was made by 
moving the probe in about 15 discrete steps, stopping at 
-each step to acquire data. Because of the probe noise men-, 
tioned earlier, at each step the signals were electronically 
integrated tor about 3 seconds. Then the integrated signals 
and all system parameters were recorded and the probe 
was moved to the next position. Because of the blow (town 
configuration, this procedure limited the number of points 
per traverse to about 15. 

The test apparatus, instruments, data accuracy, 
operating procedure, and probe noise are discussed in de- 
tail in (25). 

RESULTS 

Prior to the actual data runs adiabatic profiles in 
both gas and near-critical nitrogen were obtained upstream 
of the heater plate. These were integrated tor mass flow 
rate and the results compared very favorably with the val- 
ues from both flow meters (fig. 2). These measurements 
also indicated that the profile approaching the plate was 
flat tor a distance of 10 mm above the plate. Profiles 
measured over the plate in both adiabatic and heated nitro- 
gen gas exhibited power law behavior. The adiabatic near- 
critical profiles over the plate also exhibited power law be- 
havior. This is all documented in (25). This preliminary 
intormation was used to establish confidence in the meas- 
urements. 

In order to obtain an overall view of the results of 
the present study, it is of interest to examine the velocity 
and temperature profiles of a near-critical fluid in oom- 
pariaon to the same fluid under the same condittons except 
at free stream temperatures away from the critical tem- 
perature. These data are presented in Figs. 5 and 6. The 
thermodynamic ranges covered by these three runs are 
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those denoted on the density curve, Fig» 1. In Fig. 6 the 
low density, or **gas, ** and the high density, or *liquid»’ 
velocity profiles fall virtually on top of each other and are 
in the l/6th to lAth power range. The near-critical pro- 
file on the other hand is quite unusual In appearance, ex- 
hibiting a maximum near the wall. This is certainly con- 
sistent with the M-shaped profiles measured by Wood (4, 5) 
and Wilson (6) in pipes. The Grashof number associated 
I with the near -critical profile is substantially greater than 
the others. This suggests that the unusual velocity profile 
is the result of the large density gradient across the 
boundary layer producing a strong body force effect, even 
though the Reynolds number is close to one million. This 
is consistent with the results of the recent upflow and 
downflow tube tests (9-13). The near-critical temperature 
profile is considerably flatter than those tor the same 
fluid away from the critical point. 

After illustrating the unusual nature of near-critical 
flow fields, it then becomes of interest to examine the in- 
fluence of various system parameters on near-critical ve- 
locity and temperature profiles. 




In Figs. 7 and 8 the influence of free stream temper- 
ature on near-critical velocity' and temperature profiles is 
illustrated. At the lowest temperature, 89® K (Tq = 

126.2® K), the velocity profile appears quite conventional 
in shape. An increase of 10® K toward the critical tem- 
perature yields a profile that is virtually flat over the ma- 
jority of the boundary layer. Several of these »*flat” ve- 
locity profiles occurred in the course of the experiment 
and they will be commented on later. A continued in- 
crease of free stream temperature 10® K toward the cri- 
tical temperature yields a velocity profile with a maximum 
near the wall, similar to that displayed in Fig. 5. Fur- 
ther increases in free stream temperature toward the cri- 
tical temperature produce velocity profiles which appar- 
ently peak closer and closer to the wall. The peak cannot 
be recorded in these cases because of the limitation on 
probe travel. Corresponding to these velocity profile 
changes is a change in Grashof number of a factor of 15. 


The temperature profiles, Fig. 8, become increas- 
ingly flat as the free stream temperature approaches the 
critical temperature. This is consistent with Wood^s re- 
sults (4). On the other hand the trend in heat transfer co- 
efficients is not consistent with Wood's. He showed a 
strong increase in heat transfer coeHicient as the bulk 
fluid approached the critical temperature. The present 
results appear to indicate a decrease* The difference ap- 
pears to lie in the fact that the present results were ob- 
tained at constant Reynolds number, while Wood's results 
were at constant velocity (4). Observe that at constant 
Reynolds number the velocity changes a factor of 2 in the 
present experiment. This suggests that the question of 
whether heat transfer coefficient increases or decreases 
near the critical point is somewhat a function of the path 
along which it is examined. 


After reviewing the influence of free stream temper- 
ature on the velocity and temperature profiles the question 
arises as to whether or not the profile behavior is strictly 
a function of the free stream conditions. This can be ex- 
amined by holding all the free stream conditions - particu- 


larly temperature - constant while varying surface heat 
flux. The results of this Investigation are shown in Fig. 9. 

It can be seen from Fig. 9 that varying heat flux and holding 
free stream conditions constant produces the same velocity 
profile behavior as varying free stream temperature 
(cf., fig. 7). 

Two facts are common to both sets of data. When the 
unusual velocity profiles occur, the wall temperature is 
above and the free stream temperature is below the criti- 
cal temperature and the Grashof number has increased sig- 
nificantly. When the unusual velocity profiles occur the 
conditions are such that the sharp property changes asso- 
ciated with the critical point occur within the boundary 
layer. The condition of the wall-to-free stream tempera- 
ture spanning the critical temperature can be achieved 
either by increasing the free stream temperature towards 
the critical point holding heat flux constant or by increasing 
heat flux holding free stream conditions constant. 

The above tests were all nm under a condition of con- 
stant free stream Reynolds number. Figure 10 shows the 
effect of varying Reynolds number holding free stream 
temperature and surface heat flux constant. This effec- 
tively varies the forces in the flow field mechanically as 
opposed to thermodynamically as in the above cases. 

Again toe complete range of profile shapes can be pro- 
duced. This result is particularly significant because it 
indicates that sufficiently high inertia forces can diminish 
the influence of the strong density change. NoMce that the 
wall temperature in all cases in Fig. 10 is above the cri- 
tical temperature, making the density gradient large and 
the Grashof number quite high. Still toe Reynolds numbers 
are high enough in some of the cases to return the profiles 
to a somewhat normal shape. 

The central theme then of these Fig. 7 to 10 is that the 
veloc y profiles can be discussed in terms of conventional 
interactions of viscous, inertia, and body forces as ex- 
pressed in standard nondimensional munberS) toe Reynolds 
and Grashof numbers. The fact that the flow field respemds 
in an e?q)ected manner to variations in conventional non- 
dimensional parameters suggests that analyses based on 
turbulent mixing, such as suggested by Deissler (15), 
would be productive. The incorporation of a variable den- 
sity in the manner of (16-20) would appear to be important. 
The data indicates that toe inclusion of the body force term, 
as has been done by Hsu and Smith (17), is necessary, even 
at very high Reynolds numbers. 

It has been pointed out that the parametric relation- 
ship Gr/Re^ is considered significant in assessing the 
role of body force in near-critical heat transfer (9-11,17). 
The determination of when body forces are significant re- 
quires a certain arbitrary element in toe choice. In the 
tube experiments the demarcation was established as the 
point whore there was no significant differences in heat 
transfer between upward and downward' flow.- For the mo- 
ment, in toe present experiment, the appearance of pro- 
files which are "flat" to within 0.25 mm offer a possible 
hinge point. Five such profiles occurred in the experi- 
ment; three appear in Figs. 7, 9, and 10. They occur in 
a range: 
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1.6 <Gr^/Re^<3.0 

The actual criterion should probably be lower, since when 
the profiles are **flat»* body forces are already significant. 
However, since the ’’flat*’ profiles occur in a narrow 
range of Grx/Re|, they probably can serve as a convenient 
and significant reference condition. The spread in 
Grx/Re| above is not strictly random scatter. There ap- 
pears to be a heat flux influence; however, there is no 
clear trend. It may be, however, that . Gr^/RCx is an 
oversimplification . 

While the evidence appears strong that body forces are 
significant in the present results, other forces, such as 
axial acceleration, could also be important. Tests per- 
formed in downward or horizontal flow could be definitive 
on this question and would be a logical extension of the 
present experiment. 

In a few instances it was possible to obtain high speed 
(5000 frame/sec) movies of the near-critical boundary 
layer during a probe traverse. One set of these are shown 
along with the corresponding velocity profile data in 
Fig. 11. The view is across the plate; the flow is upward; 
and the probe is traversing right to left. Each picture is a 
single frame from a roll of film, which was taken with the 
pxobe stationary at the position indicated by the solid cir- 
^ cles. The plate surface cannot be seen since the back- 
I lighting would not penetrate the boundary layer. The ap- 
1 proximate position of the wall and the probe tip is marked 
‘ with arrows in the photographs. Note that as the probe 
t nears the wall it disappears into the boundary layer. 

The photographs indicate a highly agitated boundary 
layer. They could be interpreted as a high level of turbu- 
lence. On the other hand, the photographs could be cited 
as similar to "film boiling. " These photographic results 
I are qualitative and open to multiple interpretattons. 

^ * CONCLUSKDNS 

\ 

The main objective of this work was to make detailed 
vebcity and temperature profile measurements near the 
thermodynamic critical point. The geometry was kept 
simple both for subsequent analysis and to avoid complex 
flow interactions, thus highlighting the influence of the 
. strongly varying near-critical properties. A sufficiently 
. wide range of parameters to illustrate the various flow pat- 
terns was examined. Typical data have been presented 
^ herein. 

When compared to nitrogen in a state condition away 
the critical point but under the same operating conditions, 
the velocity profiles in near-critical nitrogen are quite 
unusual frequently exhibiting a maximum in the region be- 
tween the wall and free stream . 

Systematic investigation over a wide range of near- 
critical operating parameters reveals that velocity profiles 
over a vertical flat plate in the critical region exhibit a 
’ strong body force effect. It appears that the property vari- 
ations, particularly density, are such as to make body 
forces significant even at very high Reynolds numbers, 
i Body forces are significant where Grx/Rcx Is on the 


order of one. j 

The near-critical temperature profiles become in- ; 

creasingly flat as the free stream fluid approaches the cri- 
tical temperature. 

The fact that the vebcity profiles respond in an ex- j 
pec ted manner when traditional parameters, such Reynolds ^ 
and Grashof number, are varied suggests that the fbw ^ 

field can be characterized by turbulent fbw models. The ' 
turbulent fbw models in the literature wh*..i include den- 
sity gradient and body force effects are suggested as the « 
first avenue of analysis of the present data. 
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Figure Z - Flow system and test sedion schematic. 
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Figure 3. - Heater plate assembly. 











(b) Vlewat43‘^ 
angle from the 
tip (approx. 8 
times actual 
size). 


Overall view (actual size). 



(c) View looking directly at the 
tip (approx. 21 times actual 
size). 


Figure 4. - Actual photographs of combination pitot* 
static and thermocouple probe. 










Poo/Pc - 1. 127 q » 15. 35 watts/cm 



Photographs of the near-critical boundary layer. 










